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We have tried to measure \V u b\ using inclusive semileptonic decay measurements with three 
different methods. The full reconstruction tag is the one of the most feasible method to 
measure inclusive process. Semileptonic B — > iv tag is another method to separate signal 
and tag B mesons. New kinematic variable is used to enhance double semileptonic decays of 
both B meson. We also improve the neutrino reconstruction technique solving a combination 
of decay products into signal and tag B mesons using the simulated annealing method. 



1 Introduction 



The precise measurements of \V u b\ is one of the potential measurements to show the first dis- 
crepancy of KM-schemed standard modeffl However current |V^| measurements are limited by 
statistics, model related systematics and theoretical uncertainties. B factory experiments have 
an opportunity to overcome statistical limitation easily and the latter two uncertainties are main 
theme of the analysis. Exclusive analysis has introduced inevitable theoretical uncertainties on 
hadronic part though the analysis is simpler. Contrary inclusive analysis has difficulties due to 
ambiguous kinematics but has a chance to reduce theoretical uncertainty, if larger kinematic 
region could be measured. Theorists propose formula to extract \V u b\ from a partial branching 
fraction well defined kinematic region using the Operator Product Expansion(OPE) based on 
Heavy Quark Effective Theory(HQET) or Heavy Quark Expansion (HQE) method^. Measure- 
ments of the shape function, hadronic mass moments and lepton spectra moments may help to 
reduce theoretical uncertainties. 

From experimental point of view, the difficulties of inclusive studies are based on a lack 
of knowledge related to (l)whether event is BB or continuum, (2)which particles are decay 
daughters of either B meson, (3)whether the signal B decay into u£u or civ (4)ambiguity of 



kinematics. The neutrino reconstruction technique can solve the last issue with an assumption of 
all produced particles are observed except neutrino in the final state, though detection efficiency 
will be sacrificed. The complete full reconstruction of the opposite B meson in hadronic decay 
mode must be the best solution for the issue (1) and (2). The combination of two techniques 
is an ideal method to study B — ► X u tv and it is equivalent to the complete reconstruction of 
two B meson from T(45). However the efficiency of the full reconstruction method is not so 
large due to variety of B meson decays in hadronic modes. Belle has tried to develop other new 
methods to study inclusive b — > uiv decay besides the full-reconstruction method. 

We tried to include the semileptonic decay modes as a B meson tag as they are well estab- 
lished and has larger branching fractions as a decay mode. The semileptonic decay mode for B 
meson tag has statistical advantage against the full reconstruction method but an introduction 
of another neutrino is inevitable in the final state in addition to v from the signal b — > uiv. Using 
all constraints, kinematics can be solved with a small ambiguity as far as all decay products are 
measured except two neutrinos. Newly introduced cut parameter which define a existence of 
kinematic solution works to select BB dual semileptonic decay events. 

The neutrino reconstruction method is improved by using a constraint from B meson mo- 
mentum which is enabled by a classification of decay particles into daughters from the tag B 
meson or from the signal B meson, in order to reduce time for calculation. Combination prob- 
lem is solved minimizing probability density functions with the simulated annealing method in 
order to reduce time for calculation. As tag B meson is selected without specifying decay modes, 
reconstruction efficiency can be larger and this method has a statistical advantage against other 
methods. 

2 Full reconstruction method 

Once one of B meson is reconstructed in hadronic decay mode, events are BB from T(4S) and 
rest of all particles are assigned to the other B. Though we are not ready yet to show the result 
of b — > uiv study, this method provides almost unbiased B decay sample which is very useful 
for inclusive studies. Here we show the results on inclusive lepton study from B meson( a study 
of \V c b\ ) as a superior feature of the full reconstruction. 

B meson is reconstructed with following hadronic decay modes, 

£0 d* + it- , D* + p- , D*+a^ , D+-K- , D+ p- , D + a^ , 

B~ -> D* TT-,D*°p-,D*°a{,D 7r- 

as a tag side B meson. We did not include decay modes having larger background events as 
further selection in kinematics is not applicable in lepton inclusive study. We obtained 22k 
B° tagged events and 25k B " tagged events from 78/6 -1 on-resonance data and 8.8/6 -1 off- 
resonance data taken at Belle® detector with KEKB asymmetric collider^!. Number of events 
are estimated from the fitting of AE( Energy difference) distribution after selecting on M&( 
beam-constrained B meson mass ) distribution. Lepton signal is checked from the rest of tracks 
with momentum larger than 0.6GeV/c. Merits of this method are that the flavor and charge 
of the signal B meson is known from the tagged B meson, which enable to obtain inclusive 
branching fractions for B° and B~ meson separately and inclusive lepton spectra from B meson 
rest-frame, not from Y(4S) rest-frame. 

Preliminary results of inclusive lepton spectra for B° and B~ meson and for electron and 
muon are measured after background subtraction and corrections. Averaged lepton inclusive 
branching fraction is B(B — > X c iv) = 0.1119 rfc 0.0020 exp . ± 0.0031theo. , which is consistent with 
the previous results^. The ratio of charged B to neutral B inclusive branching fraction is 
B{B + -> X c £v)/B(B° X c iv) = 1.14 ± 0.04 exp . ± 0.01 tft , eo ., which is consistent with the life 




time ratio of B meson. This branching fraction is converted into |V^| using the formula based 
on the operator product expansion(OPE) and the heavy quark expansion(HQE)^ and result in 
\V*\ = (4.13 ± 0.07 exp . ± 0.25^*0 x 1(T 2 . 

Application of the full reconstruction to the b — ► utv is in progress and will be shown in 
near future. 

3 flWfo tag 

The second method is using D^&v tag. As these semileptonic modes are the most established 
decay modes in B meson decay and have large branching fractions, these are good for B meson 
tag. So we handle two B mesons decaying like Bt ag — > D^*'lu and B s i g — > X u £u. Separation of 
two B mesons are automatically done if B tag is specified and rest of particles except signal lepton 
in the event are assigned to X system in the signal side. But we have to treat two neutrinos in 
the final state which is non-trivial thing. 

In the single B meson which decays into the semileptonic mode like B — > Dlis, we observe D 
system and £ and also we know the absolute momentum of B meson and its mass at S-factory 
experiment on T(4J5). An assumption of zero-mass for neutrino derives the angle between 
observed D£ system and B meson's direction(^s 1 ), which means the true B meson direction is 
on the surface of the cone defined with angle 6b± around the direction of D£ system. The missing 
degree of freedom(DOF) in this kinematics is only one and corresponds to the azimuthal angle 
of B direction. When we extend this situation to two B decays from T(4S I ) and both B meson 
decay into semileptonc modes (B\ — > Dlv and B2 — > X£u), we have two missing DOF in total 
but they are filled out from the constraint that two B mesons move back-to-back each other at 
T(45) rest frame. This implies the real B direction is on the cross lines of two cones defined by 
each B semileptonic decay systems when one of B meson decay system is spatially inverted. If 
the unit vector of Di and X£ are defined as ni(0, 0, 1) and 712(0, sin #12, cos #12) (#12 shows angle 
between D£ and X£ system), real B's directional vector is tib(xb, Vb-, z b) where zb = cos0b 13 
Db = (cos 0b 2 — cos 8b 2 cos #12)/ sin 612 and 



Xfl = ±W1 ; — (COS 2 #Bi + COS 0d o — 2 COS 0Bi COS 9bo cos #12). 

V sin 0i2 

If the assumption of two B decay into semileptonic modes is correct and all decay products are 
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Figure 2: The Mx distributions( solid square for data and histogram for generic MC ) for the all-charged mode 
after all requirements are shown in the top figure. The background from B — > X c £v subtracted Mx distribution is 
shown in the bottom. Solid squares show for data and histogram shows uiv MC based on ISGW2. Right figures 

show the same for the 7r°-associated mode. 



observed except two neutrinos, x\ (inside of square root of xb) must be greater than and this 
provides severe selection condition of both B semileptonic decay events. 

3.1 Validity of Dlv tag 

The validity of this method is checked using B — > X c iv decays in signal side where tag side is 
also decaying to D^iv. This sample is obtained from 78/6 -1 data ( corresponding to 85 million 
BB events ) just requiring X system contains a single charged K in addition to existence of 
and oppositely charged lepton pairs whose momentum exceed l.QGeV/c where rest of particles 
are assigned to X system, total charge of event is zero and x 2 B > —0.2. Neutral energy is required 
to be less than 300MeV . Clear and sharp D and D* peaks are observed in the Mx distribution 
and smaller tails are very consistent between data and MC which is normalized by number of 
BB events, where MC doesn't include b — > u components. Lepton spectrum is also consistent 
between data and MC. This method can select X c tv events with high purity and good Mx 
resolution. 



3.2 B -> X u iv 

In order to select X u tv we require X system doesn't contain charged K, the other requirements 
are same as for \V c b\ sample selection. When total neutral activity in X system is less than 
300MeV, events are assigned to all-charged mode. When the neutral energy exceeds 300MeV 
and 7r°s are reconstructed from neutral clusters, reconstructed ir° is attached to X system and the 
event is assigned to 7r°-associated mode. Rest of events are discarded as poorly measured events 
to keep good Mx resolution and S/N. Fig. [21 show Mx distribution of B — > X u tv candidates 
for all-charged and ir° associated modes. Sharp peaks at D and D* are attributed from D^*' — > 
Kgnir;Kg — > ir + ir~ as we didn't remove K$. Data show clear excesses in Mx distribution 
comparing to b — > c MC distribution at tt£v in all-charged mode and at piv in both modes. If 
the signal region of Mx is defined below 1.5GeV/c 2 , we obtain 82 ± 19 excess in all-charged 
mode and 92 ± 21 in 7r°-associated mode, where statistical error includes both data and MC 
statistics. In exclusive decay mode, we observe 30 ± 6.4 events in tt~ £ + v and 28 ± 18 events 
(depend on non-resonant shape) in p°£ + u. 
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Figure 3: The efficiency corrected Mx distribution for data( soiid square with error ) and for the hybrid model( 

histogram ). 

3.3 Branching Fraction 

In order to obtain the branching fraction and V u b, corrections are very important and largely 
depend on models we choose. Corrections are divided into two parts, one is a detection efficiency 
and the other is a correction of extrapolation to unobserved kinematic regions. We use 7 different 
models to investigate two corrections. The exclusive model (IS GW2)Ld and the hybrid model 
combining exclusive model below u> in M x and inclusive model ^El above it are chosen for 
the detection efficiency. Correction factor of extrapolation is estimated by inclusive models with 
several sets of b-quark mass and ^ values. The maximum differences are assigned to systematic 
error of these corrections. 

Branching fraction is calculated by 

B = l^ N obs (M x )-N BG (M x ) I 1_ 

2 e(M x ) e ext 2N BB ' 

where N j, s and N B g show observed and b — > c MC expected number of events, e(M x ) shows a 
detection efficiency of i-th M x bin and e ex t expresses a correction for the extrapolation. Obtained 
branching fraction is 

B{B -» XJu) = (2.62 ± 0.63 sta t. ± 0.23 sys . ± 0.05 6 ^ c ± 0.41 b ^ u ) x 1(T 3 . 
4 Advanced neutrino reconstruction 

For the standard neutrino reconstruction method, a neutrino is reconstructed as a missing 
momentum in the final state and it cannot distinguish BB events from continuum events. This 
advanced method separate decay daughters into the signal side and tag side using the probability 
density function defined by various kinematic variables. There are a bunch of combinations to 
separate decay products into two categories and it is not practical to try all of combinations. 
The annealing technique^ is introduced to obtain a solution of the combination problem within 
a reasonable iteration. After classification of decay daughters, B meson constraints can be used 
to calculate kinematic variables such as the B meson momentum, B meson energy and so on, 
which improve the signal to noise ratio significantly. 

The annealing method is a algorithm to solve the combination problem. We define the 
incorrect probability for a solution using the probability density functions for a correctly sep- 
arated sample and randomly separated one by MC events where probability density functions 





Figure 4: The density probability functions for B meson momentum in CM-rest frame( p* B ), B meson energy in 
the rest frame( Eg ), the missing mass square for neutrino( MM 2 ) and the incorrect probability for the random 
sample ( before minimization ) are shown in the left and the same figures after minimization are shown in the 

right. 



are defined by several kinematic distributions. 

w _ PDF (random; p* B , E B , M M 2 , ...) 

~~ PDF(correct;p* B ,E* B ,MM 2 , ...) + PDF (random; p* B ,E* B ,MM 2 , ...) 

Search for the optimal combination is carried out by minimizing W through changing a paricle 
from signal B category to the tag side category and vice versa. With a certain probability a 
assignment of a particle is not taken even it provides better W in order to prevent a stack into 
the local minimum. Fig. |I] show several kinematic variables before minimization and after. 

After interation of this process, the optimal combination for two B meson is obtained. We 
select events whose W is consistent with zero, obtained B meson candidate is consistent with 
the real B meson using P B and E* B and missing mass calculated with B meson constraint is 
consistent with neutrino mass. 



4-1 Validity of this method 

This new method is checked using two different control samples extracted from real data and 
MC data. Mx distribution is studied using B — > D*tv sample which has been selected requiring 
a existance of D* from the mass difference of M(Dtt) and M(D) and cos9b-d*£ an d D * iv is 
considered as the signal B meson. Processing this neutrino reconstruction scheme to the sample, 
two B mesons are separated and the most probable B{— > Xtv) meson is selected which must be 
same as the B meson tagged in sample selection. Observed Mx spectrum has a peak at Md* and 
the shape is consistent between data and MC except for its normalization( e(data) / 'e(MC) = 
0.902±0.040 ). q 2 distribution is also checked using B — > J/^X sample where one of leptons from 
J/* decay is treated as neutrino. Obtained q 2 distribution shows a peak at q 2 = M 2 ^ and same 
shape for data and MC, though the normalization between data and MC is slightly different like 
in the Mx distribution from D * Iv samples ( e(data) j 'e(MC) = 0.862 ± 0.046 ). Normalization 
differences between data and MC are considered to be from uncertainty of hadronic decay mode 
in B meson decay table of MC generation and this effect (e(data)/e(MC) = 0.886 ± 0.034 ) is 
corrected in the X u £u analysis. 
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Figure 5: Mx vs. q 2 distribution for selected sample( a ), the projected Mx distribution after requirement of 
q 2 < 7GeV 2 for data( cross ) and MC( histogram ) ( b ) and the background subtracted Mx distribution ( c ) 
are shown. In the last figure, two solid histgrams show expected distributions from inclusive models with two 
different mj ( mj = 4.95GeV/c 2 and mj, = 4.Q5GeV/c 2 ) and shaded histograms show expected distributions from 

ISGW2 model for -nlv and piv. 



4.2 Branching fraction of B — > X u lv 

Through this advanced neutrino reconstruction, B — > Xlv candidate are selected from 85 Million 
BB sample corresponding to 78/6 -1 . Selected sample shown in fig. G3-(a) contains huge X c tv 
components which are considered as backgrounds for B — > X u iv study. The kinematic regions for 
the signal are limtted in Mx < l.SGeV/c 2 and q 2 > 7GeV 2 /c 2 in order to reduce X c bacgrounds 
and to optimize the signal sensitivity. We obtaine 1148 events excess to the background in these 
region. The branching fraction of B — > X u iv is obtained from 

1 (N obs -ax N bkg ) 1 

D — > J\. u tV — — 



2 axe// 2N BS 

where N obs and N bkg is obtained number of events in the signal region for data and MC respec- 
tively a = e(data)/e(MC) and e// shows a detection efficiency of signal estimated from MC. 
Models used for efficiency estimations are same as the previous method. Obtained branching 
fraction is 

B(B -» X u lv) = (1.64 ± 0.1W ± 0.36 SJ/S . ± 0.28 6 ^ c ± 0.22 b ^ u ) x 1CT 3 . 
5 Extraction of \V ub \ and systematics 

We used the following formula^ to extract from B(B — » X u £u), 

\V ub \ = 0.00445 x ( ^ -, u£u)1.55ps )1/2 ^ ± ± Q 



where r b = 1.606 ± 0.012 is used as a average of B + and B° lifetime^. Extracted \V ub \ are very 
prepliminary and \V ub \ = (5.00 ± 0.60 siai . ± 0.23 sys . ± 0.05 6 ^ c ± 0.39 6 _ u ± 0M theo .) x 10~ 3 ( for 



Table 1: Systematics of B(B — » X u iv) for D'*'£v tag and v reconstruction method. 
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£>*^z/ ), \V ub \ = (3.96 ± 0.17 siai . ± 0AA sys . ± 0.34 fe ^ c ± 0.26 6 __ n ± 0.29 theo .) x 10~ 3 ( for advanced 
v recon. ). Obtained values are consistent each other within error as well as results from other 
measurements. Systematic errors from several items are considered and listed in tabled for each 
method. 

6 Conclusion 

We present preliminary results on B — > X u lv branching fraction using two new methods. 

- for D*iv 

B(B X u lv) = (2.62 ± 0.63 siat . ± 0.23 sys . ± 0.05 6 ^ c ± 0.41 6 _ u ) x HT 3 , 
\V ub \ = (5.00 ± 0m stat . ± 0.23 sys . ± 0.05 fe _ c ± 0.39 6 _ u ± o.m theo .) x 10~ 3 

- for advanced v recon. 

B(B -> X u lv) = (1.64 ± 0.14 stQi . ± 0.36 sys . ± 0.28 6 ^ c ± 0.22 6 _ u ) x 10~ 3 , 
\V ub \ = (3-96 ± 0.17^. ± 0.44 sys . ± 0.34 6 _ c ± 0.26 6 ^ u ± 0.29 theo .) x lO" 3 

Bothe methods provide good S/N in addition to good Mx resolution for D^lv tag and good 
efficiency for the advanced neutrino reconstruction tag. 
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